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The influence of two temperature 
generalized thermoelastic diffusion inside a 

spherical shell 
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Abstract- In this investigation, the problem of 
elasto-thermo-diffusion interaction insider a spherical shell is 
concerned in the context of the mechanism of two temperature 
generalized thermoelasticity theory. The inner and outer 
boundaries of the spherical shell are traction free and 
subjected to heating. The chemical potential is also assumed to 
be a function of time on the boundary of the shell. The problem 
is based on the theory of two temperature generalized 
thermoelastic diffusion with one relaxation time (i.e. two 
temperature Lord-Shulman (2TLS) model). To obtain the 
general solution the integral transform technique is used. The 
solution obtained in the Laplace transform domain by using a 
direct approach. The inversion of the transformed solution is 
carried out by applying the method of Bellman. The numerical 
estimates for the thermophysical medium are obtained for 
copper material and have been shown graphically. The 
influence of diffusion on thermoelastic stresses, conductive 
temperature, thermodynamic temperature, displacement, 
concentration, chemical potential inside the shell are observed 
for Lord-Shulman model. The results in the absence of 
diffusion are also found as a particular case. 

Index terms - elasto-thermo-diffusion interactions, generalized 
thermoelasticity, thermoelastic diffusion, two temperature 


I. INTRODUCTION 

Diffusion can be defined as the migration of the particles 
from region of high concentration to region of lower 
concentration. The recent interest in the study of this 
phenomenon is due to its extensive applications in 
geophysics and many industrial applications. In most of 
these applications, the concentration is calculated using 
what is known as Fick’s law. This is a simple law that does 
not take into consideration the mutual interaction between 
the introduced substance and the medium into which it is 
introduced or the effect of the temperature on this 
interaction. 

Biot [1] develop the coupled theory of thermo-elasticity to 
deal with defeat of the uncoupled theory that mechanical 
cause has no effect on the temperature field. In this theory, 
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the heat equation has a parabolic form which predicts an 
infinite speed for the propagation of mechanical wave. 

The theory of generalized thermoelasticity with one 
relaxation time was introduced by Lord and Shulman [2]. 
This theory was extended by Dhaliwal and Sherief [3]. In the 
theory, the Maxwell-Cattaneo law of heat conduction 
replaces the conventional Fourier’s law. For this theory, 
Ignaczak [4] studied the uniqueness of solution; Sherief [5] 
proved uniqueness and stability. The fundamental solution 
was developed by Sherief for spherical regions and by 
Sherief and Anwar [6] for cylindrical region. 

Nowcaki [7] - [10] developed the theory of thermoelastic 
diffusion. In this theory, classical coupled thermoelastic 
model is used. Later on, Gawinecki et al [1 1] proved a theory 
on uniqueness and regularity of the solution for a nonlinear 
parabolic thermoelastic diffusion problem. For the same 
problem a theorem about the global existence of the solution 
was established by Gawinecki and Szymaniec [12]. Sherief 
et al. [22] and, later on, Kumar and Kansal [23] introduced 
the generalized theories of thermoelastic diffusion in the 
frame of LS and GL theories by introducing thermal 
relaxation time parameters and diffusion relaxation time 
parameters into the governing equations, which allow the 
finite speeds of propagation of waves inside the medium. 
Sherief and Salah investigated the problem of a 
thermoelastic half space in the context of the theory of 
generalized thermoelastic diffusion with one relaxation time. 
Singh [26], [27] studied the reflection phenomena of waves 
from the free surface of an elastic solid under theory of 
generalized thermodiffusion. Aouadi [28]-[33] also gave 
some attention on thermoelastic diffusion and generalized 
thermoelastic diffusion. Othman et al. [33] reported some 
studied on the effects of diffusion on a two-dimensional 
problem of generalized thermoelasticity in the context of the 
mechanism of Green-Naghdi theory. The theory introduced 
by Sherief et al. [22], Kothari and Mukhopadhyay [38], 
presented the Galerkin-type representation of solutions for 
thermoelastic diffusion. In the context of the same theory, 
variational and reciprocity theorems have been established 
by Kumar et al. [38]. 

The linearized version of the two-temperature theory (2TT) 
has been studied by many authors. Chan and Gurtin [13] and 
Chen et al. [14], [15] have formulated a theory of heat 
conduction in deformable bodies, which depends on two 
distinct temperatures (a) the conductive temperature ([) and (b) 
the thermodynamic temperature 0. Lesan [17] has 
established uniqueness and reciprocity theorems for 2TT. 
The existence, structural stability and spatial behavior of the 
solution in 2TT have been discussed by Quintanilla [18]. 
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The key element that sets the two-temperature 
thermoelasticity (2TT) apart from the classical theory of 
thermoelasticity (CTE) is the material parameter % (>0), 
called the temperature discrepancy [15]. Specifically, if x=0, 
then <|>=e and the field equations of the 2TT reduce to those 
of CTE. It should be pointed out that 2TT suffer from 
so-called paradox of heat conduction, i.e., the prediction that 
a thermal disturbance at some point in a body is felt instantly, 
but unequally, through-out the body. 

The present paper is concerned with the investigation of 
disturbances in a homogenous, isotropic temperature 
dependent elastic medium with two temperature generalized 
thermodiffusion. The analytical expressions for the 
displacement components, thermoelastic stresses, 
conductive temperature, thermoelastic temperature, 
concentration and chemical potential are obtained in the 
physical domain whose boundaries are traction free and 
subjected to a time dependent temperature and chemical 
potential in the context of 2TLS model. The Laplace 
transform technique is used to obtain the general solution. 
To get the solution in the physical domain, the inversion of 
the transformed solution is carried out by applying the 
method of Bellman. The influences inside the shell is 
analyzed for a copper like material and depicted graphically. 
The effect of two temperature parameter in the presence and 
absence of diffusion are analyzed theoretically and 
numerically. The most significant points arising out from 
our analysis are highlighted. 


II. MATHEMATICAL FORMULATION OF THE 
PROBLEM: 


We consider an isotropic homogenous thermoelastic 
spherical shell with inner radius a and outer radius b in a 
uniform temperature T 0 . Let the body be referred to 

spherical polar coordinate (r, 6 , cp) with the origin at the 
center O of the cavity. Since we consider thermoelastic 
interactions which are spherically symmetric, so all the 
functions considered will be function of the radial distance r 
and the time t only. It follows that the displacement vector 
hi thermodynamic temperature 0 and Conductive 
temperature (f) have the following forms: 

u = (u(r,t), 0,0), 6 = 6{r,t), (p = (p{rf (1.1) 


In the context of two temperature generalized thermoelastic 
diffusion based on LS theory, the equation of motion, the 
equation of heat conduction and the equation of mass 
diffusion in absence of body forces for a linearly isotropic 
generalized thermoelastic solid are, respectively, 


8 2 u & „ 86 8C 

P^2=( A + 1 V)^-P^-P2^ 

ot or or or 

KW 2 </> = pc E (6 + t 0 6) + T 0 {3 1 (e + r 0 e) + 

cT 0 (C + t 0 C) 


( 1 . 2 ) 

(1.3) 


Dp 2 V 2 e + DcV 2 0 + C + t°C = DdV 2 C, 

where 0 is the thermodynamic temperature, (j) is the 


conductive temperature, X and p are Lame's constants, 


p is the density, C is the mass concentration, /?, , /? 2 are 
the material constants given by 

/? 2 = (3X + 2p)a c , = (3X + 2p)a t , 

in which a t and a c are, respectively, the coefficient of 

linear thermal expansion and linear diffusion expansion, K 
is the thermal conductivity; D is the diffusion coefficient and 
C E is the specific heat at constant strain, r 0 is the thermal 

relaxation time andr 0 is the diffusion relaxation time, c 
and d are the measures of the thermo-diffusion effect and 
diffusive effect, respectively. V 2 is the Laplacian, given by 

v>=LU| r ^ 


r 2 dr 


dr 


The strain components are given by 
du u 

e rr ~ T - ’ e ee ~ e H ~ ~ ’ 

dr , r 

and thus the cubical dilation will be 

du u 1 d 2 
e = — + 2 — = — — (ru ) . 
dr r r dr 

The constitutive equations are given by 
de 

a rr =2p- + Xe-p i e-p 1 C 

or 


(1.5) 


( 1 . 6 ) 


(1.7) 


( 1 . 8 ) 

Where 


&ee =^h> = ?-Li- + ^e- P x O- P 2 C 
r 

P = -fi 2 e + dC-cO , (1.9) 

P is the chemical potential per unit mass of the diffusive 
material in the elastic body, cr. is the stress tensor. 

The relation between the conductive temperature (j) and the 
thermodynamic temperature 6 is given by, 




( 1 . 10 ) 


Where X (>0) is the two temperature parameter. 

For convenience, the following dimensionless quantities are 
used: 

u = cpju, r - c^pr, 9 f = 

Pi c jj- M 


P x 0 


C' = 


</>'= 

(A + 2 p) (A+2 n)‘ 

[rV;] = Ci 2 77[r°,r 0 ]. 


p'=A 

A 


(A+ip) 

t' = c 2 rjt, 


2 (A + 2 u) pc F 

Where c x = and Tj ~ 


P 


K 


Then, the governing equations are given by equations (1.2) - 
(1.4) and equations (1.7)-(1.10) can be expressed in the 
following forms (dropping the primes for convenience): 


(l.H) 


d 2 u _ de 56 8C 
dt 2 dr dr dr ’ 

VV = (6 + t 0 $) + s(e + T^e) + scc x (C + TqC\ (1.12) 


(f> — 6 = aN 2 (j) 

V 2 e + a 1 V 2 0+ a 2 (C + r°C) = a 3 V 2 C, 


(1.13) 

(1.14) 
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4 u Q 
2 2 u 

Cng = (1 T )e + ---0-C, 

ee pi ' p r 


P = —e + a 3 C - a,6. 


Where Cl, = 


cpc. 


PiPl 

T„P “ o2 


a , = 


pc. 


V D0 : 


a, = 


(1.15) 


(1.16) 
(1.17) 
d pci 


2 ’ 3 Pi 


_ u. ?2 o2 _ ^ + 

~ ,2„2 ’ A' ~~ 


P C \ C E 


P 


G> = %cW- 


0 = 


1 v sa,coa 3 - scoa 




C 


1 + coci 3 1 + 1 + coci 3 

where <2 3 = v(l + r 0 s). 

Applying divergence operator on (1.21), we get 

(V 2 -/)e = V 2 C + V 2 0. 


Using the equation (1.22), equations (1.12)-(1.14) and 
(1.23) becomes 


V 2 (f> — lyfi + LjC + L^c , 


(1.24) 


[(1 + !,«») V 2 e ] - s 2 e = L V 2 </> + (1 ■ - Lp)' V 2 C, ( 1 .25) 

d 3 

(1 - L.cojV 1 ! + ^V 2 0= [(a 3 + coa ,1^ )V 2 - 


a 2 s(l + r°s)]C, 
where, 


A 


a 3 


1 + 


(1.26) 


r a 3 a x s r a 3 s 

^2 ~ — ’ “3 ~ — 5 

1 + ci 3 CO 1 + Cl^CO 


III. BOUNDARY CONDITIONS: 

From the physical phenomena for the problem, we will take 
the boundary conditions as: 

(1) The boundaries of the shell are assumed to be traction 
free, i.e.: 

cr rr (r,0 = 0 On r = a,b fort >0. (1.18) 

(2) Boundaries of the shell are subjected to a thermal shock 
in the form 

(f) = (f) l H(t) on r — a, t> 0, 

= </> 2 H(t) on r — b, t > 0. (1.19) 

(3) The chemical potential is also assumed to be a known 
function of time at the boundaries of the shell, that is: 

P = P l H(t) on r = a, t > 0, 

= P 2 H(t) on r = b, t >0. (1.20) 


IV. SOLUTIONS IN LAPLACE TRANSFORM 
DOMAIN: 

Applying the Laplace transform defined by the 

oo 

relation, f(r,s) = | f(r,t)e~ s, dt Re(s) > 0, 

0 

to equation (1.11) and using homogeneous initial conditions, 
we get 

2 - de d6 dC 

su = . ( 1 . 21 ) 

dr dr dr 

Now applying the Laplace transform on equation (1.12) and 
using (1.13), we get 


^=e-y~-e-c, 

p r 

2 2m-- 

P = -e +a 3 C-a,0. 

Now, using Equations (1.24)-(1.26), we obtain 
( V 6 - by + by - b 3 ) (e , 4>, C) = 0. 

Where we used the notations 


d-27) 

d-28) 

d-29) 

(1.30) 


b,= 


1 alil-l^cof-aliapI^apXl +L l sco) 
2 „. /, , o 


[a,L l a 3 L 3 {l-L 1 o))+ 


oz,L,u 3 Lj2(1 + + ci 3 $2^(1 1 - f >t)(l I - LfSco) + L,a 3 ( 1 + L,sco) 

( a 3 + Lflp) + Lfi] (1 + L,£(o){a 3 + Lp,®) + a]s 2 ( a 3 + + 
^L i a 3 (a i +bap)+^I^a 3 (\-L 2 (o)-L^al(\-L 2 co) ^ ] 


(1.31) 


b 2 = — t — t — [a 3 L, 

^(l-Z^ta) -a 3 (a 3 + L 2 a l co)(l + L l £co) 

2 2 2 2 0 

L,a 3 s + L,a 3 s (a 3 + L 2 a,6)) + a 3 sa 2 (l + T s)L, (1.32) 

(1 + L,£(o) + a 2 3 s 2 a 2 (1 + rV) + a 2 L,L 3 a 3 s{\ + r°s)] 


[dap (\ + t ] s)L s \ 


3 2 / 


a l 3 (1 - Pp) f - a 2 (a 3 + L 2 a 1 to)( 1 + L,£co) 


(1.33) 


Equation (1.29) can be factorized as 
(V 2 - kl )(V 2 -k 2 2 )(V 2 -k 2 3 )(e,0,C) = 0, (1.34) 

where k x ,k 2 , k 3 are the roots with positive real part of the 
characteristic equation 

k 6 -b,k 4 +b 2 k 2 -b 3 =0, 


( 1 . 22 ) anc j are gj ven by [ 24 ] 


(1.23) 
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+ 4 

■1 [2p sin q + b x , 


f [A - cos q 1 + sin t/)] 


!~ [A + p(V3 cos q - sin qr)] , 

where 



(1.35) 

(1.36) 
d-37) 


P=$~\ Q = 


sin 1 v 


v=- 


26 ( - 96 ^+ 276 , 


2/ 


(1.38) 


Therefore, the solution of equation (1.30), which is bounded 
at infinity, is given by 

j(r,s) = ^^[A l (s)I v (kf)+B t (s)K v (k l r)], (1.39) 

sir £l 72 /2 

e(r,s) = ^2^[A!(s)f 1/ (k i r)+B l '(s)K v (k i r)], (1.40) 

sir tT /2 /2 

C(r,s) = ^[A;W v (kf)+B;{s)K v (k i r)l (1.41) 

s/r “T /2 /2 


where 4 , 4' , A! ", B n B[, B" (for i = 1, 2, 3) are 
parameters depending only on s. Iy 2 is the modified Bessel 
function of the first kind of order 1/2 and Ky 2 is the 

modified Bessel function of second kind of order 1/ 2 . Now, 
from Equations (1.24)-(1.26) and (1.39)-(1.40), we obtain 

the following relations: A r = — A , A" = — A , B r = — B , 

^ l i l 7 l j l J l 1 l 

di d t a i 

f. 

B"=^B, (1.42) 

* d. 1 

where, 


p t = a 3 (1 - L^scd)^ + \LJ^ - (1 - L 2 coj]kf , 

d i = k 2 [a 3 Z^ (1 - L^co) + a 3 Z^ (1 + sl^co)] - s 2 a 3 L 2 , 
f { - a 3 (1 + sbfiS)^ - [a 3 Z 1 (1 + sl^co) + 5 , 2 a 3 

+Z^Z^]/^ + s u 3 L^. 

Thus we have 

^r,s) = ^ g f[A i (s)I v (k i r)+B i (s)K v (k i r)], 

Vr m A /2 /2 

C(M) = -j=E%4W M+B^K M]. 

sir £1 A /2 /2 


(1.43) 

(1.44) 


Where 

^3 0 L 2 sco)k i + [Z^Z^ ^ 3^4 (1 Z*a>)]*, ? 

^ 2 i = d 3 (l-\- sL^co)ki — \cl 3 L ^ (1 + sL^cd) + s a 3 + L^L^Ik- + £ a 3 Ly. 
Using the relation between u and e from (1.5) and (1 .43), we 


get the solution for the dimensionless form of displacement 
as follows: 


u{r, S ) = ^j]^[A i ( S )I y (k i r)-B i ( S )K y (k i r)]. (1.45) 

Vr m Mi 72 72 


Therefore, from Equations ( 1 ,27)-( 1 .29), (1.39) and 
(1.43)-(1.45), we get 

3 . , . _ _ ... *+ U: _ 

IT 

(1-46) 


^Ar,s) = ^j]A i (s)[gJ v (k i r)-^^I y (k i r)]- 

slrtt /2 P rkpi /2 


(1-47) 


(r, s) = 4= £ A (AkJ y (k,r)+ 2 ^' I y (*,r )]+ 

VrM /2 P rk i d i n 

is+s^K (k,r)--^-K (*,r)], 

“7 72 P rk i d i 72 

P(r,s) = ^j^g 5 XMs)I v (k i r)+Bps)K l/ (k i r)]. (1.48) 

sir tr 72 72 


Where 
f = £l_i 

? 4 4 ’ 

Pi aJ, 

C, ^ 

5 ‘ d, d. 




To evaluate the unknown parameters, we shall use the 
Laplace transformation of the boundary condition 
(1 . 1 8)-(l .20), together with equations (1.39), (1.46) and 
(1.48); we have the following set of six linear equations in 
six 

unknowns: 


EA {s)[g v I yikfl)- 


*Pi 


1=1 

3 


(3 2 ak i d i 7i 


- 1 3/ (kj a )] + 




+ - 


4 Pi 


(1.49) 


1=1 


pak i d i 7i 


K 3/ (k.a)\ = 0 , 




2>,wte MP) 


+ - 


4 Pi 


i = 1 


Fbk.d. 7i 


Ky(kb)\={), 


X [A 0)/ 1 / (M) +B, (s)K v (k,a)\ = 


1=1 


g[A(^ W+^,- W1 = 
E^[A(M 1/ (A«)+AW^i/(M)] = 


i=l 


fasfa 

? 

</> 2 sfb 

Pisfa 

5 

P 2 s[b 


(1.50) 

(1.51) 

(1.52) 

(1.53) 


X? 5 ,[A(M 1/ /(M)+AW^(M)] = ^^- d-54) 


1=1 
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We can obtain A^y), A 2 (y), A 3 (y), B x (s), B 2 (s ), 
B 3 (s) by solving the above linear system of equations 

(1.49)-(1.54). This completes the solution of the present 
problem in the Laplace transform domain. 


IV. SPECIAL CASE: (WITHOUT DIFFUSION) 


By putting C — 0 and j5 — 0 into equation (1.2), (1.3), 
(1.7) and (1.8) and neglect the diffusion effect by eliminating 
Equations (1.4) and (1.9), we get the equations for 
conductive temperature, displacement and the stresses 
without the effect of diffusion. In this case, after some simple 
computations, equation (1.30) reduces to, 


Where a, - 


(V 4 -a 1 V 2 +a 2 )(e,(/>) = 0 . 

ci^Ly (1 + LySco) + s cl 3 + 
a 3 (l + L^sco) 
Equation (1.55) can 
as(V 2 -^)(V 2 -fc 2 2 )(^) = 0, 


(1.55) 


, a 2 


L^s a 3 


a 3 (l + 1 ^ 8 ( 0 ) 
be factorized 
(1.56) 


where k x , k 2 are the roots with positive real part of the 

characteristic equation k A — a x k 2 +a 2 = 0 . 

Using the solution of equation (1.55), given by 

0(r,s) = ^[A i (s)I v (k i r)+B i (s)K y (k i r)], (1.57) 

Vr i~i 7i /2 

e(r,s) = ^j][A:(s)I v (k i r)+B;(s)K v (k i r)], 

7 r i=1 22 22 (1.58) 

where A j , A- , B j , B, (for i =1, 2) are parameters 

depending only on s, we get 

& 


i 2 


Vr nM, /2 /2 

2 


4ft 


(1.59) 


<*n ( r > s ) = -7= Z A' Mb'* V*i r ) - -tfTjh/ (V)l + 

Vr m /2 P 44 h 


Y^mi^vAkr)- 


4ft 


prk i d i h 


KM l 




(1.60) 


^B i (s)[g,,K y (k / )--^K y (k i r)}. 

m /2 4 44 h 


(1.61) 


Where 


A 


Pi = LJc. . d i = a 3 [(l + I^sco)^ -s ], $- 3 , =-^-1, 

d ; 


2l)Pl 
P 2 d, 


$4i = (l __ 3r) J -1- 


In view of the boundary conditions (1.18) and (1.19) A i (y) 
and 5(v) involved in Equations (1.57), (1.59)-(1.60) can 
be determined by the solution of following 

Z A 1/ (*,«) - A-/y (*,«)] + 


equations: *=i 


/3 ak i d i /2 

4 Pi 


(3 ak i d i /2 


W 3/ (fta)] = 0, 


«[&*,/ W+ JPtd K Yi <W1 = °’ 

(f> x 


1=1 

3 


XtA WA/(M+A A)^i/(M] = 


(=i 

2 


(1.63) 


(1-64) 


X[4(X>A/ (k^+B^sWy (kjb)} = 
i=1 /2 /2 J (1.65) 

This completes the solution of the present problem in the 
Laplace transform domain. 


V. NUMERICAL RESULTS AND DISCUSSION: 

In order to illustrate theoretical results in the preceding 
sections, we now present some numerical results. To get the 
solutions for the displacement, radial stress, shear stress, 
conductive temperature, thermodynamic temperature, 
chemical potential and mass concentration in the physical 
domain, we have to apply Laplace inversion formula to the 
equations (1.49)-(1.54) respectively. Here we adopt the 
method of Bellman et al. [45] for inversion and choose a time 
span given by seven values of time t { , i =1 to 7 at which u , 
<r rr , cr^, P and C are evaluated from the negative of 

logarithms of the roots of the shifted Legendre polynomial of 
degree 7. For the illustration we consider copper material 
with material constants. The physical data in SI units for 
which given as follows [24] : 

2 = 7.76xl0 10 yV.»f 2 , // = 3.86xl0 10 (V.m" 2 , p = m4kg.nf\ 

K = 386W.m- 1 K~\ c E = 383.1 J.kg' l K A , T 0 = 293K, 

a c =1.98xl0“ 4 m 3 kg~\ a, =1.78xl0' 5 K“ 1 , 
s = 0.0168, c = 1.2x10VkV 2 , 

d = 0.9 x 10 6 m 5 kg A s" 1 , a x =5.43, a 2 =0.533, a 3 =36.24. 

Also we have taken </> 0 =1, T (j = 0.01, r°=0.1 and R = 1 
for computational purposes. 

Figures 1-6 representing the variation of displacement, 
radial stress, shear stress, conductive temperature, chemical 
potential and mass concentration along the radius of the 
sphere for both the one-temperature Lord-Shulman (1TLS) 
(for 0 = 0 )and the two-temperature Lord-Shulman (2TLS) 
(for co = 0. 1 ) theories. The computation were carried out for 
both large time ( £=0.35 ) and small time (t =0.026). The 
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variation of the field is observed when the step input of 
conductive temperatures with (j) x —\ and (j) 2 = 1 applied on 
the inner boundary a=l and outer boundary b=2 respectively. 
The step input of chemical potential with P x = 1 and P 2 = 1 

are also applied on the boundaries of the shell. Our 
computations were also carried out in absence of diffusion, 
by using Equations (1.57), (1.59)-(1.61). In figure 1-6 
dashed lines represent the case in the thermoelastic diffusive 
medium in the context of 1TLS theory ( 0 = 0 ) (1TLS WD). 
Doted lines represent the distribution in thermoelastic 
diffusive medium in the context of 2TLS ( co = 0. 1 ) (2TLS 
WD) theory. The solid lines and thin lines correspond to the 
thermoelastic distributions (that is, in the absence of 
diffusion) under the 1TLS theory (0 = 0) (1TLS WOD) and 
the 2TLS theory ($ = 0.1) (2TLS WOD) respectively. 



Figure 1 : Distribution of u r at t=0.026 (a) and t=0.35 (b). 

Figures 1(a) and 1(b), depicts the variation of the 
displacement component (u) against radial distance (r) 
inside the spherical shell for both small time ( t =0.026) and 
large time ( t =0.35), respectively. From figure 1(a) it is 
observed that, a mild effect of diffusion is in the 
interval 1 < r < 1.2 . Moreover, the displacement component is 


compressive in nature near both the boundaries at lower 
time. The graphs of displacement under with diffusion and 
without diffusion for 1TLS (0 = 0) and 2TLS ( 0 = 0.1) 
theories are almost merged. As may seen from the figure 
1(b), it is observed that the radial displacement (u) takes 
negative values near the inner boundary ( 1 < r < 1 .7 ) of the 
sphere for both the cases with diffusion (WD) and without 
diffusion (WOD) for both 1TLS model and 2TLS model. The 
absolute value of this field also increases with the increase of 
time in both the media. 



1TLSWD 

2TLSWD 

1TLSW0D 

2TLSW0D 


(a) 



1TLS WD 
2TLS WD 
1TLS WOD 
2TLS WOD 


(b) 

Figure 2: Distribution of G rr at t= 0.026 (a) and t=0.35 

(b). 

Figures 2(a) and 2(b) are plotted to show the variation of 
radial stress ( <y rr ) against radial distance (r) inside the 

spherical shell. At both boundaries the radial stress is noted 
to be zero, which also agrees with the theoretical boundary 
conditions. At time t =0.026, the radial stress CT rr is 

compressive in nature throughout the medium (1 < r < 2 ), 
whereas with the advancement of time it becomes fully 
compressive at the middle zone for both theories. For both t 
=0.026 and t =0.35, the influence of diffusion is significant 
towards the middle of the shell under the 1TLS (0 = 0) and 
the 2TLS ( co = 0. 1 ) theories. 
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Figures 3(a) and 3(b), the variation of shear stress ( cr ^ ) 

against radial distance (r) inside the spherical shell are 
observed. At both small and large times, it may seen from the 
figure, the hoop stress is fully compressive in nature in all 
cases. At small time the effect of diffusion is very prominent 
inside the shell, whereas, with the increase of the time it 
decreases under both 1TLS (<y = 0) and 2TLS (^ = 0.1) 
theories. 




(b) 


Figure 4: Distribution of (j) at t=0.026 (a) and t= 0.35 (b). 

Figures 4(a) and 4(b) display the variation of temperature (j) 
against radial distance (r) inside the spherical shell. We can 
see that the magnitude of temperature field shows the 
maximum value at both the boundaries. At small time, with 
the increase of radial distance towards the middle it 
decreases and becomes minimum at the middle of the shell. 
For both larger and small time the difference of numerical 
value for both with diffusion and without diffusion under 
1TLS ( co = 0 ) and 2TLS ( co = 0. 1 ) theories are also noted. 




Figure 5: Distribution of P at t=0.026 (a) and t= 0.35 (b). 


Figures 5(a) and 5(b) represent the variation of chemical 
potential ( P ) against radial distance (r) inside the spherical 
shell for the thermoelastic diffusive medium. At both 
boundaries the chemical potential is noted to be in 
agreement with the boundary conditions. The difference 
between the chemical potential are more prominent under 
the 1TLS ( co = 0 ) theory as compared to the 2TLS ( 0 = 0. 1 ) 
theory. 
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(a) 



Figure 6: Distribution of C at t=0.026 (a) and t= 0.35 (b). 

Figures 6(a) and 6(b) shows the variation of concentration 
(C) against radial distance (r) inside the spherical shell for 
the thermoelastic diffusive medium. Like temperature field 
the mass concentration shows the maximum value at both 
the boundaries and becoming minimum towards the middle. 
The 1TLS (<y = 0) theory predicts a significantly different 
trend as compared to 2TLS ( 0 = 0. 1 ) theory. 


VI. CONCLUSIONS: 

The problem of investigating the thermoelastic 
displacements, stresses, conductive temperature, 
thermoelastic temperature produced in a homogenous 
isotropic spherical shell of the thermoelastic diffusive 
medium is stated in the context of 1TLS and 2TLS theories. 
We also compare our results with the corresponding results 
in thermoelastic medium (without diffusion). According the 
analysis above and the numerical results presented in 
figs. 1-6, we can conclude the following important points: 

(i) The presence of diffusion plays an important role in all 
the quantities. The influence of diffusion is more significant 
on displacement, radial stress and shearing stress, as 
compared to conductive temperature and thermoelastic 
temperature fields. At small time, the influence of diffusion 
on stresses near the boundary is more significant. With the 
increase of time the region of influence shifts towards the 
middle. 

(ii) The significant difference is also noted in the physical 
quantities for one temperature and two temperature LS 
models. Two temperature theory is more realistic then the 
one temperature theory in the case of both generalized 
thermoelasticity with diffusion and without diffusion. This 
study is very important for microscale problems, because in 
these cases the material parameters are temperature 
dependent. 
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